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ON THE KERNEL FUNCTION OF THE INTEGRAL EQUATION RELATING THE LIFT AND 
DOWNWASH DISTRIBUTIONS OF OSCULATING FINITE WINGS 

IN SUBSONIC FLOW ’ 

By Charles E. Watkins, Harry L. Bunyan, and Donald S. Woolston 


SUMMARY 

This report treats (he kernel function of an integral equation 
that relates a known or prescribed downwash distribution to an 
unknown lifi distribution for a harmonically osciUating finite 
wing in compressible subsonic flow. The kernel function is 
reduced to a form that can be accurately evaluated by sepayating 
the kernel function into two parts: a part in which the singular- 
ities are isolated and analytically expressed and a nonsingular 
part ivhich may be tabulated. The form of the kernel function 
for the sonic case {Mach number of 1) is treated separately. In 
addition j results for the special cases of iMach number of 0 
{incompressible case) and frequency of 0 {steady case) are given. 

The derivation of the integral equation which involves this 
kernel function, originally performed elsewhere {see, for example, 
A'. id Technical Memorandum 979), is reproduced as an 
appendix. Another appendix gives the reduction of the form of 
the kernel function obtained herein for the three-dimensional 
case to a knoivn result of Possio for two-dimensional flow. A 
third appendix contains some remarks on the evaluation of the 
kernel function, and a fourth appendix presents an alternate 
form of expression for the kernel f unction. 

INTRODUCTION 

The analytical determination of air forces on oscillating 
wings in subsonic flow lias been a continuing problem for the 
past 30 ATars. Throughout the first and greater part of 
this time, efforts were directed mainl 3 ^ toward the determina- 
tion of forces on Agings in incompressible flow. These efforts 
iia\e led to important closed-form solutions for rigid wings 
in tAA’o-dimensional floAA" (ref. 1), to solutions in terms of 
series of Legendre functions for distorting AAungs of circular 
plan form (refs. 2 and 3), and to man}" approximate, yet 
useful, results for AAungs of elliptic, rectangular, and tri- 
angular plan form (see, for example, refs. 4 to 12). 

Although these results for incomprossiblo flow play a 
highly significant role in applications of unsteady aerod}"- 
namic theory, the ad\"ent of higher and higher speed aircraft 
during the last 15 }"ears has brought a groAving need for 
knowledge of the effect that the compressibilit}' of air might 
have on unstead}" air forces, or for analytically deriA’ed un- 
steady air forces based on a compressible medium. The 
transition to results for a compressible fluid from those for 
an incompressible fluid is not likely to be accomplislied by 


applications of simple transformations or correction factors, 
such as the well-knoAvn Fraud tl-Glauort factor for steady 
floAv. This difficulty is associated Avitb the fact that tlio time 
required for signals arising at one point in the medium to 
reach other points gi\"cs rise not only to cliaiiges in magni- 
tudes of forces but also to additional phase lags between 
instantaneous positions, A"elocities, and accelerations of the 
AA"ing and the corresponding instantaneous forces associat(*d 
AAUth these quantities. In order to obtain results for the 
compressible case, it therefore appears necessary to deal 
directly with the boundar}"-A"aliie problem for this case. 

The boundary-A'alue problem for a tAA"o-dimonsional Aving 
in compressible floAV has been successfully attacked from two 
points of A"icw. First, by consideration of an acceleration or 
pressure potential, Possio (ref. 13) reduced the problem to that 
of an integral equation relating a prescribed downwash dis- 
tribution to an unknoAAm lift distribution. The kernel of tliis 
integral equation, AAdiicIi is a rather abstruse function, Avas 
reduced to a form that, except at singular points, could be 
evaluated. Scliwarz (ref. 14) later isolated and determined 
the analytic heliavior of the singular points of Possio’s results 
and made fairly oxtensiA^e tables of the kernel function. 
These tabular values AA'ore used by various investigators 
(for examples, refs. 15 and IG) to obtain, by numerical 
procedures, initial tables of force and moment coefTicients 
for oscillating AA'ings in compressible subsonic flow. 

T)ie second successful approach to the solution of the 
bouTulary-A’alue problem for a two-dimensional wing (see 
refs. 17 to 19) is acliieA'ed by a transformation to ellijitic 
cooi’dinates followed by a separation of A'ariables tliat reduces 
tile boundary-A"aluc problem from one in partial-differential 
equations to one in ordinary differential equations of the 
Matliieu type. Tlie solutions turn out as infinite series in 
terms of Matbieu functions. Xiimcrical results obtained 
recently by this procedure agi’ee with results jireAuously ob- 
tained by tlie numerical procedures using tlie kernel func- 
tion (see, for example, ref. 20). 

With regard to boundary -A-ahie problems for finile Avings 
in compressible flow, it appears that the procedure of sepa- 
ration of A"ariablcs could be a feasible approacli only for 
AA’ings of A’cry special plan forms siuh as a circle or an ellipse. 
In any ease, tlic dcA’elopmeiit of the appropriate mathe- 
matical functions for a particular plan form would become 
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lii^^hly iiivolvoil. Oji the other lunul, it appear? that approxi- 
mate procedure? similar to those used for two-dimensional 
wing:? might afford an approach to solution? of these prob- 
lems wliicli, thougli laborious, might be liandled by routine 
numerical methods. 

The kernel function of the integral equation relating pres- 
sure and downwash for tlie thive-dimensional case appears 
as an improper integral. The purpose of tliis report is to 
treat and discuss this kernel function. The improper integral 
is reduced to a form that can be aeeui‘at(‘ly evaluated by 
numerical procedures. The form and order of all its singular- 
ities are detei'mimal and an expression for the kernel function 
is derived in ^^■hieh th<‘ singularities are isolated. Special 
forms of the kernel for the sonic case {^f= Ij, the incompres- 
sible case (d/==0), and the steady case (^'-=0) are presented. 

series expansion in powers of the reduced-frequency param- 
eter k is developed. 

The availability of the kernel in a form which can be 
rapidly eA'aluated makes possible the use of numerical pro- 
cedures, similar to those used in tlie two-dimensional case, 
to obtain aerodynamic forces for finite wings. 


ANALYSIS 

TNTEGR.\L EQUATION AND ORIGIN.AL FORM OF KERNEL FUNCTION 

The main purpose of this analysis is to ti'cat the kernel 
function of an integral eipiation that relates a known or 
prescribed downwash distribution to an unknown lift dis- 
tribution for a barmonically oscillating finite' wing in com- 
pressible subsonic flow. The inti'gral equation referred to 
can be obtained by emplo^'ing the Prandtl acceleration 
potential to treat linearized boundary-value problems foj* 
oscillating finite wings by means of doublet distributions. 
De rivation of this integral equation from tlie linearizeul 
boundary-value proldem for a wing is a preliminary taslc 
that has been done elsewhere (see, for example, ref. 21), hut 
it is reprodu(“ed heroin as an appendix for the sake of com- 
pleteness. 

In keeping with the concepts of linear theory, tlie wing is 
considered a plane impenetrnhle surface S which lies nearl,y 
in the j’?/-])hmo as indicated in sketch 1 : 
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K{j'o,yo) 

K'{xo,Vo) 
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velocity of sound 

ITankel functions of second kind of zero 
and fii*st order, respectively 
modified Bessel functions of first kind of 
zero and first order, respectivel3" 

Bessel function of first kind of zero order 
modified Bessel functions of second kind of 
zero and first order, respectively 
kernel function of integral equation 
singular part of K{xo,yo) 
redueed-froquencA’ parameter, loo/V 
modified Struve functions of zero anti first 
order, respect h^el 3’ 
unknown lift distribution 
reference length 
^fach number, T7c 
pressure 

region of 3-?/-plane occupied b\' wing 
time 

forward velocit.y of wing 

amplitude function of prescribed downwasb, 

Cartesian coordinates 


Euler^s constant 

velocit3’ potential 
acceleration potential 
fluid density 

circular frequenc}’ of oscillation 



Skoldi 1. 

The x,y,z coordinate S3’stem and the surface S" are assumetl to 
move in the negative j-diretUion at a uniform velocit3' V. 

In terms of these coordinates, tlie integral equation maA' be 
formall3' written as 


w (x, y) Jj L{^,ri) K (xo, 2/0) dk dv ( 1 ) 

iS 

where w{x,y) is the amplitude function of the prescribed 
downwasb, .K(xn,ya) f, y— 17) is the kernel function 

and pln^sicalh^ represents the contribution to downwash at 
a field point {x,y) due to a pulsating pressure doublet of unit 
strength located at any point (J,?;), and L ({,??) is the unknown 
lift distribution or local doublet strengtli. 
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The kernel fiinelion may he malhenial ieally dermed by the 
following improper integral expression (sec eq. (A12), 
appendix A) : 


52 Cr 




(/X (2) 


where M is Alach number, /3=^'l— d/% w=co/F/3‘, oj is the 
circular fi’cqiiency of oscillation, T" is the velocity, and 
X is the variable of integration. Evaluation of this in logical 
constitutes a main difTiculty in obtaining aerodjmamic 
coefficients for oscillating finite wings in compressible flow. 
The present analysis is therefore devoted to reducing it to a 
form that can be accuratcl 3 ' evaluated by numerical pro- 
cedures combined with the use of tables of certain tabulated 
functions. The form and order of all its singularities arc 
determined, and an expression for the kernel function is 
derived in which the singularities are isolated. 

REDUCTION OF THE KERNEL FUNCTION 

In considering the reduction of the kernel function 
K{xo,yo), the integral involved can, for convenience, bo 
written as the sum of two integrals, namel.y 


j: 


(, (x — X/ X2-j-r2j 


<l\^- 
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{) X^+r2) 


/X"b 


d\ 


Therefore, 


K{j„,y„)=\im v ;2«~ ’’ (F)=lim . (F.+Fj) 

:-*[) z-^0 


where 


and 




— fu) (^XF-A/ X^Fr^) 


(/X 


F2 = 
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Jy ^Jw(x-.U^ X2 Ft% 




(I\ 


(3) 


(4) 




( 6 ) 


and where r= 

The integrals Fi and Fo are treated separately’ in succeeding 
sections. The final forms arc given in equations (15) and 
(19), respectively’. 

Evaluation of Fi. — The integral Fi can be converted to a 
form that can be more easily' handled by" writing 


'■=x; 


p — iZ) -A / Vx^F r' 
g-i3X ' .■■■ ■■ ■■ ■ — (l\ 


and introducing the folloAA'ing relation (see p. 41G of ref. 22) 


, = MT\) _ TdT 

vV+r" Jo 

= Jo(TX) " T dT- 


f' -ij ill p — (r r* 

-X 


(7) 


In the first integral of tliese last two integrals, make the 
substitution 


and in (he second integral make the substitution 




Then 


j- = I ' Jo(X-\V^ i A/-W-’) </r— 

^X“+r Jo 


j; 


f"'' JD(X’v’-iF’:;=-r-) ih 


(s^ 


(It is of interest to note, in the expression on the left of eq, 
(8), that X and r appear in the same manner. Tlio roles of 
these two quantities could, therefore, be int err hanged in 
the expression on the right.) 

With use of equation (8), the equation for Fi can be 
Avritten as 

ii’,= ,ZX f-" Jo(Xa't*+- 1/-:3-) (/r- 

; e-" J„(X ^ AP5W) <ZtJ (9) 

Changing the order of integration in each integral (Avhicli is 
a legitimate step because the integrands involved satisfy 
the continuity conditions required for such operations) leads 
to the following expression for F\\ 

F,= J” C-" dr I^J^” f-“‘ Jo (X rfxj- 

[Jo" * f/x] ( 10 ) 

The integrals within the brackets in equation (10) may be 
evaluated from tables of Fourier or Laplace transforms as 
(see, for example, pair no. 55 of appendix III of ref. 23) 

f“ Jo (X^F+M^) d\= ,-A^ 

Jo \T — p‘W“ 

f" Jo (Xv-l/^l?^) d\= 

Jo r -j-p cj 

SO that 

f” -J==dr^ 

Jo Jo \T^-T^~C0' 

The first integral in equation (11) can be written as 

f" r dr-- i 


( 11 ) 


or 


r ^S^Z= dr= f 

J« Jo 


eosh $ JQ. 

0 2 , 


iX” ‘ 
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The first integral on (he right of equation (11a) is given on 
page ISl of reference 22 as 


i: 




where /lo is the modified Bessel function of tlie second kind 
of zero order. The second integral on (lie right of equation 
(11a) is given on page 33S of reference 22 as 

where lo is tlie modified Bessel function of tlie first kind of 
zero order and Lo is the modified Struve function of zero 
order. Then, the first integral of equation (11) can be 
VTitten as 

X? 

• dr=Ko(^ccr)—^ [7o(/3wr)— Xo(/3w/')] 

[/„ (l?, v+? 

-to (y Vj/o*+2^)J (12) 

Note that the end result indicated in equation (12) is in- 
dependent of Mach number. Tlie second integral in equa- 
tion (11) may be written in anotlier form as 


g-*rT 

do ’'“Jo 


-i VffoHz-)r 




(13) 


This integral has not been reduced to closed form; however, 
it is nonsingular and can be readily handled by numerical 
methods. 


(5? L„=£fr"' 


Combining equations (12) and (13) gives the following 
expression for Fi : 

F,=iu (f [a (f 

\ V<P+2''^r 




/1 + r^ 


(H) 


By performing the difTereiitiations indicated in equation (4), 
there is obtained for the first part of equation (4) the follow- 
ing expression : 




+ 

( 15 ) 


All terms of this expression other than the integral nmy 
be evaluated at small intervals of ijo from existing tables, 
except at 2 / 0 “fi where the function is singular. The integral 
is well behaved and can be accurately evaluated hy numerical 
or approximate procedures. The type and order of the 
singularities at yo=0 are discussed in a later section. 

Evaluation of — In order to reduce the integral F 2 , 

equation (6), it is convenient to make the substitution 


so that 


F2 = 


/: 


\ — r sinh d 

. -lio 


(sInh 0—3/ cosh 0) ^ 


( 10 ) 

( 17 ) 


Noting that z appears only in r and performing the differen- 
tiations indicated in equation (4) ^fields 


(sinh 0—3/ cosh 0) 


Xo 


yo\'xo^+P'yo‘ 


fiiw (jo-3/ 




Xo 


, I w (id - 3/ 


J *si 

0 


fi\Vo] 


[0^ cosh 0 — (cosh d~M sinh 0)] ei^^[yoI(slnh0-3/cosh0) 




[ ^ I w ( 7 * - 3/ Vro-'-;?* Ve’) — ^-i«3//?lvo|] 1^-P— r cosh 6 (sinh 0-3/ cosh 0 ) 

M[l/o|Jo 


(18) 


or, by reverting completely to Cartesian coordinates througli 
equation (16), there is obtained 



This expression vanisbes, as it should, for Xy=0 and, like 
that in equation (15), has singularities at ?/o=0 wliieh, also, 
will be handled in a later section. The integral that remains, 


like the integral remaining in equation (15), is nonsingular 
and simple in form and can be readily evaluated by numerical 
procedures. 

Expression for the kernel in terms of nondimensional 
length variables. — Equations (15) and (19) can now be 
combined to gi^^c a reduced form of the kernel function 
^(xo,2/o)« Howci^or, in application, the 3"ariables Xq and y^ 
are emplo 3 Td, for convenience, in nondimensional form. 
This is accomplished h^^ considering these variables in a 
new sense to mean tlmt thc}^ have been referred to some 
chosen length I and by introducing the rcduccd-frequencj’ 
parameter k—IcojV, The variables will be used in this new 
sense tlnoughoiit the remainder of the report. The kernel 
can be vTilten in terms of these nondimensional variables as 
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iMk]yo\ r.\r/0 

“ -Jo ^ 


\(kxof+^\kyof fi.Jo-.'\/V(to)=+«=ui/o)=] 
i p^'o ^[x-.vVv+^(tM>] 

^W2AO“Jo 


+ 


fZxj- 


( 20 ) 


An alternate and perhaps more desirable form of expression 
for the kernel function is given in appendix D. 

Note that this expression for K(iQ,y^) can be considered as 
a funetioi) of only three parameters, namely, k\yo\, kx^, and M. 
To be more specific, the first two terms are functions only of 
^‘ll/ol; the next two terms are functions of and d/; and 
the last two terms are functions of Jclyal, and M. 

Equation (20) constitutes the principal residt of this report. 
Some partial checks as to its correctness arc: (1) For A*=0, 
it reduces, as discussed subsequently, to tlu* downwash of a 
pressure doublet in steady flow and (2) an integration Avith 
regard to the y-direction between the limits — «> to + <» 
yields Possio's result for the two-dimensional case. This 
integration is carried out in appendix B. Other special forms 
of the kernel function for A/=l, A/=0, and k=0 are derived 
in subsequent sections. A power series expansion of the 
kernel which is applicable for certain ranges of the parameters 
^'jyol, kxQ. and M is presented. In the section immediately 
following, the orders and types of tlie singularities of the 
kernel function are discussed. 


DISCUSSION OF THE SINGULARITIES OF THE KERNEL FUNCTION 

As previously indicated, the kernel function becomes singu- 
lar or indeterminate at yo=0. The forms that the kernel 
function takes when it becomes singular are of particular 
importance in applications to lifting surface theory. It is 
therefore desirable to extract and treat the singularities 
separately. 

This extraction can be conveniently made by considering 
the value of K(xo,yo)y equation (20), at points on the semi- 
circumference of a small ellipse (sec sketch 2), the polar 
equation of which may be written as 

Xq~ c sin S 

ijq=^ cos d 




After substituting these expressions for Xq and ?yo into 
equation (20), the results may be written as 

IT/ ® / A’€ cos 0 /k€COSd\ 

9)= /v-eos^e I ^ A, j- 

iirkt cos 6 |”^ cos ^ ^ cos ^ 

ikMt coa 9 ikf (bIq 9—\f) 

e P e , ikt cose 

^ M ' 


sin 6e 


Fe * COS ^ e 


0^ 


ik C“"“ cos2 e) 


iK f 

mJo 


X ' 


/ikt cos 


\1 + T^€ ^ P ^ dr -\- 


( 22 ) 


With the use of the following series expressions for Ki(z) and 
[Ii(z)—Li(z)] (which can be obtained from ref. 22 — for A", 
see p. 80; for A, see p. 77; and for Lu see p. 329) : 


3 \4 


+ 5 ^ 4 . S-' 4 _ 
^64^1152^ 


■) 


(23) 


where y is Euler’s constant (7=0.5772157), and 


[/.(3)-i.(.)]=HJ+^-g+4+ . . . (24) 


where, because of the S 3 'mmetry of A(jo,2/o) 'vith respect to 
1/or only the limits — 7 t/2^0^x/ 2 need be examined. Note 
that in these equations values of 0 in the range — 7 t/ 2 ^ 0 <CO 
correspond to field points ahead of or upstream from the 
doublet position and values of 0 in the range 0<5^7r/2, to 
field points behind or downstream from the doublet position. 
In particular, 6=t12 corresponds to points dircctl^^ behind or 
in the wake of the doublet. 


it is found that for vanishingly small values of t the limiting 
value of the expression for K(e,d) in equation (22) is for 
A /<1 


K(e,d) 


ala 6 r _^2 ’J. p ^-,(1 -sin d) 

P \ €\l -sin ey € 2 2(1 ~M) 

~ (M-sin 0-~)]+O(e") 


(25) 
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where 0(t”) represents terms of order e” for d S 1 . Expressed 
in terms of Tu and )/o, equation (25) becomes 


K(To,yo)==—jf 

F- 


' T — (-ru +A'V+i 

I ?/oSV+/3'; 


I 

'"l/o^ ^r/+|3-!h■ 


2 “ 2(1 -M) 2 L 2 




(2G) 


Examinaliou of equation (25) shows tliat tlio ke rnel fnnet ion 
K{tfi) has singularities with respect to €= V^o*~r/3-?/i? as 
follows : 


JM- '1- [y,(0)+log e] 


(27) 


where, from equation (25), 

^ _^-(l +sin B) 


JM- 


1 —sin d 


cos^ Q 


k cos^ 6 


. 1 k(\~9>\ne) , 

/,>(0)-lo£r ^2(riM) “ 2(1— d/)(l -i-sin 0) 


(28) 


Although of no particular significance in a])plications, it is 
of interest to note that the quantities fi and /2 each ImA'e 

minimum values ^i/i|min=^ R'nd l/^Utn^log j 

— 7 t/ 2 , which corresponds to points directly ahead of the 
doublet position; and, as 9 increases from — tt/ 2 to +7 t/ 2, the 
values of these quantities continuously increase from these 
minimum values to infinite quantities as follows: 


^'(Or'is 


/=(j)l-iii» 

\-/l 5-0 
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/3"|^H-Kin (n-^) 
cos^(|-s) 
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6— U 


log 
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k cos^ 



2(1 -il/) 

1 H- sin ^ 

i-‘)] 


A‘6‘ 


(29) 


Thus A^(jo),Vo) Is singular for 6=Trj2 even when the distance 
e from the douhlet is not necessarily of zero order. This 
implies that the doublet produces a wake of discontinuous 
downwash that extends downstream from the doublet 
position to infinity. 

With knowledge of the singularities involved in the kernel 
function K{xoJJo)f an expression can be written in which the 
kernel is separated into a singular part and a nonsingular 
part (as was done bj' Schwarz, ref. 14, for the two-dimen- 
sional case) as follows 

AT/o,?/o) ^ [ AT.ro.?/o) - K' (xo4h)] + K' {x,,y,) (30) 

where K{x^,]h) is defined in equation (20) or (22) and 


K>(, ,, wC'" ^ _ 

L To'-i-jS-i/o- 

F a-o— dA'r„^+/3-)// F, ^rn^+Fyo^—^n)] 

Off 2“’=' 2(1 -d/) J 


Of ill terms of e and 6, inlrodiuT'd by oqiiations (21), 

f-*< 5 in«r ,32 ji- ^.2 

^ p L“?(rare)+T- 

Ir . /i’€(l — sin 0)1 

2 ’^s-sdzrr/) J 


( 31 ) 


(32) 


Tlie term [A(xu,?/n) — A'(>o,2/n)l in equation (30) is a continuous 
function for all values of 1% and and for values of 1/ in 
the range of O^d/^1. The term A'(ro,yo) is discontinuous 
at the douhlet position (ro = 0, ?/o=0) and at all points in 
the wake (r(,>0,?/o=0). It is to he noted, however, that 
each term of K'(xo,yo) possesses a simple indefinite integral 
witli respect to yu or with respect to ‘q=y~y^, a fact tliat 
may be usdul in some numerical applications. Tlie manner 
in wliich these integi-als arc to be evaluated is indicated in 
a siibseqiKMit section that deals with steady flow. Tlie 
limiting values at yo=0 of [A(ro,yo) — A'(ro,2/o)] for holli 
subsonic and sonic flow arc given in appendix C togeth(‘r 
with some remarks on evaluation of the kernel function. 

TREATMENT OF THE SONIC CASE 

Because of its special nature, the borderline case, ^f=^, 
between subsonic and supersonic flow deserves and requires 
separate treatment. 

As the expression for the kernel function given in 

equation (20) becomes indeterminate. It is possible, bow- 
cA^cr, to obtain conditional limiting values for the keimel by 
considering the integral F, equation (4), and breaking it into 
two integrals, Fi and F,, as was done for the general case. 

With regard to Fi, its limiting Arabic and the Auilue of its 
derivatiA^es with respect to z at 2=0 can be shown to be zero 
as d/— >1, Fiom the form of Fj giA^en by equation (14), 

£o(pvyo 

= (|, Vr')-^ [/„ (|a W+ 2=)- 




(f T“b 


But since (sec ref. 22, p. 172) 


-r“fxv^-=-A’-«(f) 

jfl VH-A 


(33) 


( 34 ) 
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and (see ref, 22, p. 332) 

;r“®'y^c/r=^[/o(.o--£o(f)i 

it may be roncluded from equation (33j that 


lim = lini 



(3o) 


(3G) 


Tile total eont ri[)ulion lo K{i\,ifQ) at therefore, arises 

from the limit of /y, equation (0), as The liinitinjr 

foiTu of Fo may he written in terms of noiulimensional co- 
ordinates as 


lim F 2 = lim 
A /— 1 ‘ 



a''X“+/3‘(?/o^+ z^) 


d\ 


(37) 


in approaelun^ the limit M=\ (from the subsonic side) in 
equation (37), it is convenient to rc])la(‘e d/by 


M=\-€ 


where e is infinitesimally small so that 


^~(l~d/) (l+d/j-€(2-0«2e 


With this approximation, equation (37) may be wi'ilteii as 


lim F 2 = lim 


r 


X — jx’(l — e) I 1-j- 






d\ 


-X 


(for j'o>0) 


(38) 


From physical (*onsiderations, the right side of equation (38) 
is to be eonsidei’cd zero for i’o = G- Th\i^ is in keeping with 
results that would be obtained if the limit under consider- 
ation were sought from theory of supersonic flow, d/]>l. 


The integral in equation (38) cannot be completely 
expressed in terms of known functions. Furthermore, since 
it is singular at its lower limit, further treatment is required 
to reduce it to a form such that its (lorivatives with respect 
to r can be numericalh’ evaluated. For this purpose the 
integral may be written as two integrals, namc'ly 


where 


and 


{F2)m.1 = F./^FF 




k--hz'i€^ 


d\ 


-=x: 




d\ 


ff/0-+2‘ 


(39) 

(40) 


(41) 


The liinils of integration in eqniition (40) are so chosen that 
the integral in this equation can he reduced to a known h)rni 
by making tlio substitution 

X=V rM- ( 2 / 0 *+ 2^— r or 


Thus, 


f" . 

Jo A 


+(yo^+^) 


H' 




,l + r= 


dr (42) 


Equation (42) may be written in terms of the integrals 
involved in F| (see eqs. (34) and (35)), namely, 


F,' =Ko (k , 2/oH --O 1 Uk A 2/u'+ :^-Uk , ?V+ r)] 

(43) 

DilTerentiating this result twice with respect to z and then 
setting 0 gives 






(44) 


DilTerentiating equation (41) twice with respect to 
setting 2=0 gives 



2 and 


(45) 


After perfoTTuing an integration by parts and collecting 
terms, equation (45) may be written as 





+ 



(40) 


Equations (44) and (4G) are combined to give 



Then, in accordance with equation (4), there is obtained 
for K{xo,ya)in.\-- 


For ro>0, 

^(•'o,yo)j/-i='^ < -XTifA'iJ/oi) 9y('T)^ 


Uk 






Vu k^y<i‘ 


+ 


A f‘"' 

k~yo'Jl;'^ 




</xj- 


and, for Xq^O, 


A(/ojyo).u=i — 0 


(47a) 


(47 b) 


The integral appearing in equation (47a) is finite and 
proper and can lie evaluated by numerical proe(‘dures. 

TREATMENT OF THE STEADY AND INCOMPRESSIBLE CASES 

Tt is of interest to eonsidei' tlie form of tin* kernel function 
given in equation (20) for some particular values of M and l\ 
In the following sections a discussion is given for the steady 
case (/' — 0) and the incompressible ease (M—0). Tln^ two- 
dimensional case is bandied in appendix B. 


391619— 5G 2 
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Reduction of the kernel for the case of steady flow, — In 
onicr (o ol)tain the reduction of tlie kernel for the case of 
steady flow, ronsid(*r the expanded form given hy equation 
(26). As ^0, there results the following expression 

A'(A,.?/o)a.o=-k( A+ , d ( 48 ) 

whicli represents the downwasli of a pressure doublet for 
steady flow. This result serves as a partial check as to the 
con*ectness of the expression for iv(/o,.Vo) given by equation 
(20). 

By replacing ?/o in equation (48) by y — ij and iiil (‘grating 
from —1 to 1 with respect to ??, there is obtained 

j_K{,,,y,)dr,=—iY J 


lim 

O 




\(k\yo\)~Li(k\y^+i 


^ , (53) 

By setting in equation (53), a form is obt aim’d which 

can be shown to agi’ee with results dei’ived by KUssner for 
the case d/=0, .Co=0 (nd. (26)). 

A SERIES EXPANSION WTTH RESPECT TO "A” 

An approximation for tlie function 

[AA/o,2/o)-A'A-ro,?/o)] 

for small values of k can l>e obtained by making use of the 
series expansions for Ki (eq. (23)) and for (Ii—Li) (eq. (24)) 
and expanding all other terms of ^(^* 0 , 7 / 0 ) (eq. (20)) into a 
power series in term.s o( k. After pimforming these expan- 
sions and collect iiig term.s with respect to powers of k, there 
is obtained for d/<^l 


where the symbol -C indi(‘ates that a principal value or A(xo,2/o)^ 

finite part of the improper inti’gral must be taken. (See, 
for (example, ref. 24 for a discussion of finite parts of sueh 
integrals.) Tliis result corresponds to the down wash pro- 
duced hy a simple horseshoe vortex two units wide. An 
equivalent expression for incompressihle flow is given, for 
example, in reference 25, where in contrast to the present 
notation, Xq has been chosen as positive forward. 

Reduction of the kernel for In ord(‘!“ to eflect I lie 

reduction of the kernel for the incompressible case, the 
expressions for Fi, equation (15), and Fo, equation (18), will 
be examined for the limit 4/— >0: 

From equation (15) 


f-ton f i3-(-vXoH;8^Vo'+Jo) , ikfi' 


e-»^ f 

\ yoS 


I H I *U/ I j 

-xVT/Sy 


2 L 


, / ] \ ?V/3- 

{7-2)-^- 




*,j,[(12J//J^-20ili^;8'+15/8«- 1 2^)y^-d2MW+ 
4 (3i\/'+6M=-l)a-„H 12|3= iM*+2AP-l)xoy^ 


and from equation (18) 


12 ^ 2 / 0 ^ log 


vxo^+isy 


"ik (* si oh ' 

lim o=r~. I smh ^ dS- 

02“ ]yo[Jo 


yoW+2/o“ 


iM^+ 10M-- l)Xo\'V+|9='?/„^- 


Integrating by parts yields 


lim f Vl/o^+X-’ e>c^d\ ^ ^ 

2/0 yo yo Jo iA) A'-^o +Vo" 

s-»0 


Combining the results from Fi and F^ gives for the kernel 
function 

A (3*0, 7 /oXw - 0= p ^ “ l~j K I (^’1 9|^ [^1 (^*1 ^0 i ) ~ 

j/oS'j-oHyo’ y« 




— 12AF/3^Xo?/o^— 




o^*(3AP— l)yoW+;8’j/o' 


For valiums of tlic paranuUc’vs that satisfy tlic following 
inequalities 


~ ( Jo- A/ A >n^-f iSW) < 1 


equation (54) yiidds results that are corn’d to within about 
2 percent. 
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Correspondingly fov M—\, equation (47) can be expanded 
to obtain 


^-ikxoj 2e 




, , To _ Tq“_£? ;/o“« ^ 

2 2;/? 2"^ 2T„^ 


1 I ^'i.Vo' , 1 To 

jLl-V-log i-+log|jj- 

I •'?To tq=' :iy»^ ypV ^ 1 

6 L 2 4yu 4 j-o 4T(,“ J 


Z Z^o 


Szo'+eyo^ log 17,”,+^" 


?/ol To' 


ipro’l 

>“ yo^e ^ I 

1= “2t„‘ J' 


To’ 3 to1/o' yo*_ yo‘ , 2/o’eJ^ 1(56) 

2 ■^Jo 12 jo-'^20toM 


For values of the pnrainelors that satisfy the following 
inequalit}': 

A-To-^2^"<2 (57) 

equation (5G) yields results that are eorre<‘l to witliin about 
2 ptTcenl. 

CONCLUDING REMARKS 

The main purpose of this report was to present the kernel 
function of the integral equation relating the downwasli to 
the lift distribution in a form tliat can be computed. This 
purpose has been achieved by the presentation of tlie kernel 
in a form given in equation (20). Tliis equation lias been 
converted to a form more suitable for calculation by isolating 
the singularities as shown in equations (30) and (31). The 
special case of M=1 is given in equations (47). The forms 
of the kernel function for other limiting cases, namely A* = 0 
and J\/=0, are given in equations (48) and (53), respectively. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 18, 1958. 
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APPENDIX A 


DERTVATIOX OF THE INTEGRAL EQUATION THAT RELATES THE DOWNWASH AND LIFT FOR A FINITE WING BASED ON 

REFERENCE 21 


In koopinj; tlu* conropts of linonr tlu*ory, tlio wing is 
considered ns a nearly ])lane impenetrable surface, this 

surface >S lie nearly in the -ry-]>lane, ns indi(*ali‘d in skcUch 1 
of the body of tlie report, and let it and the x, y, z coordinate 
system to wbicli it is referred be assumed to move at a 
uniform speed V in the negative /-direction. At tbe same 
time, let eacli point of the wing be assumed to undergo 
harmonic translations of small amplitude Zmixjj.t) at 
circular frequency co and let c represent velocity of sound in 
the medium. 

Tlie problem for an oscillating wing consists in solving the 
wave equation subject to certain boundary conditions. The 
wave equation in rectangular coordinates is 


dV ■ dV dV 



(Al) 


The independent variable ^ in equation (Al) is regarded 
herein as an acceleration potential; as such it is directly 
proportional to a perturbation pressure field and is related 
to a velocity potential 0 as follows: 


In order to complete the boundary-valm‘ problem for the 

c)0 

wing, it is desirable to calculate the downwash ic{x,y,Zff)=^^ 

associated with 0. Assuming this downwash to be harmonic 
with regard to time implies that both potentials 0 and 0 arc 
harmonic with regard to time and can be written, therefore, as 


0(r.y.r) 

^(x,y,z) 


(A3) 


With these expressions for 0 and 0, equation (A2) becomes 
independent of time and reduces to an ordinary equation 
with one independent variable, namely 

(A4) 

Tilis equation can be integrated with res])ect to x to give 
—— . ^ 

(Ao) 

where the lower limit of integration is chosen, for later 
convenience, so as to satisfy the condition that 0 vanish 
as 00 . 

10 


The ])oimdnry-\ ahu‘ problem for the wing may now be 
expressed mathematically as follows: Under the assum])tioii 
of liarmonic motion tlie differential equation, erpiation (Al), 
becomes 


5~0 ^~0 d ‘0 


C‘ 




(AG) 


In order to insure tangential flow at tbe wing surface, the 
potential must satisfy the dovmwash condition 

w(i',J/)=(5^)^_^=(' (A7) 

where w and are amplitudes of velocity and displacements, 
respectively, and are assumed to he kiiovm from the motion 
of the wing. At r = the pressure 


p=~p{4'')z^o (A8) 

must be zero at all points (/,?/) ofl’ the wing. At all points 
on the wing 0 is allowtul to be discontinuous and the value 
of p at a given point is determin(‘d by the magnitude of the 
diseontinuity in 0 at the point. In the neighborhood of the 
trailing edge, p must go to zero, corresponding to the Kutta 
condition. 

One other condition, that 0 vauisli far ahead of the wing, 
is inherently satisfied by the relation between <f> and 0 given 
in equation (A5). 

The ])otentiai 0t, al point (/, y, z) due to a harmonically 
pulsating doulilet located in the j7/-plane at (^, r\, 0) that 
salisfu^s equation (AG) is 


where 


0O = -1 


C> C L 


efi- 


] 


7?'-= vU-«f+/3^’ (;/- vr+H^z- 


(A9) 


and the factor .1 is a strength an<I dimensionality factor that 
makes possible' diff(*rent uses and interpretations of the 
])olential 0(,. If 0o is considered as an acceleration potential 
and substituted into equation (A.")), there is obtained a 
corresponding velocity potential 0o which may be Avritten as 



J 


- -J, -rfx 


(A 10) 


iiriin 
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wliorc 




The downwnsli 


9<i>o 

~6z 


associated witli ipv may be written as 


0<i>o 4 3’ 

d7~^‘^d?^ J_„ VW 


d\ 


(All) 


where XQ=x—^f w=a>/T"/3^, and With the 

use of this equation and the concept of solving linear 
boundary-value problems by means of superposition of 
elementary solutions to the governing differential equation, 
the boundary-value problem under discussion can be writ ten 
as an integral equation, namel}’ 


w (x,y)=lim 
' 2-^0 




iiilTo 



(/J dv 




(x-.VVV+rJ) 

y)?+?“ 


d\ 

(A12) 


where S represents the surface of the wing and n) 

represents an unknown lift distiibution or doublet strength 
on S. Equation (A12) may be seen to correspond essentially 
to equations (1) and (2). 

If the distribution function Z(J, r)) in equation (Al2) is 
deteiTuined in accordance with the boundaiy conditions 
discussed in the preceding paragi*aph, equation (A 12) can 
be consid(‘red as a complete solution to the boundary-value 
problem for an oscillating finite wing in compressible flow. 
It is also to be noted that equation (A12) can be considered 
to represent a solution to the so-called “indirect” problem, 
that is, that of finding the do^^mwash distribution associated 
with a given lift distribution. 


APPENDIX B 


REDUCTION OF THE KERNEL FUNCTION FOR THREE-DIMENSIONAL FLOW TO THAT FOR TWO-DIMENSIONAL FLOW 


The purpose of this appendix is to show that integration of 
the kernel function K(;XQ,yo) from — oo to + with respect 
to ij— 1 /— ?/olcads to a known result for two-dimensional flow. 
The kern(d is first modified to a form that, for the present 
case, is easier to handle. Then, after performing an integra- 
tion l)y parts on the modified kernel, the form of the kernel 
for the two-dimensional case is given (eq. (B18)). In addi- 
tion, the special cases of M=\ (cq. (B23)) and A/=0 (eq. 
(B30)) are also shown. 

The integration under consideration with respect to rj is 
equivalent to an integi'ation with respect to t/o, namely 

ij K (xo,y—v) dv=l j K (xo,yo) (hjo (Bl) 

It is remarked in advance that since 3 has been made zero 
in the expression for /iC(iro>2/o), equation (20), it is necessary 
to employ the concept of “finite parts of infinite integrals” 
when integrating this function across the singularities at 
Use of this concept gives the same results that 
could be obtained by the more arduous task of performing 
the integrations before setting z equal to zero. 

Modification of the kernel. — In order to effect the desired 
modification of the expression for K{xQ,yo) given by equation 
(20), consider the first integral of the expression, namely 

dr (B2j 

This integral can be written as 

lim— r aT+V ^/r-f F f dr 

6-^ Jo Ja///3 

(BH) 

but according to page 331 of reference 22 

Y,{d4-iJc\yo\)] (B4) 


where IIi is the unmodified Struve function of fii*st order and 
Yi is the Bessel function of the second kind of first order. In 
the limit as 5^0 these expressions have the following values: 
For the first expression in the bracket (see ref. 22, p. 329) 

lim //i(5 + fA"|yol)-//i(?‘^l2/ol) = -i:.(^*ll/ol) (Bo) 
6-^) 

and for the second expression (sec ref. 22, pp. 77 and 78) 


lim ri(6+?I*|?/ol)== — + 

i-*o 2/ 

=^^^i(^l2/ol)-/i(^il/ol) (B6) 


where IIi denotes the Hankel function of the first kind of 
first order. With the use of equations (B3) to (B6), expres- 
sion (B2) can be 'written as 




1 -{- r a'1 e Jr-|- 

jM/ff 

^ [/.a-lj/oD-iia-ii/oi)] j- (b 


7) 


Substituting this result into equation (20) of the text gives 
the modified form of KiXQ^yo) sought, namely 




ik - 






- p P 4- - f 0 — 

fi\yo\ ^My<? 




r \- 

Ja//^ 


l + r^ Jr+ 


ik p* jj(x-A/Vxj+w) 


ik p 

Myo" Jo 


d\ 


m 


Integration of modified kernel. — Since the expression for 
Kixojjo) is symmetrical with respect to 2 / 0 , that is, .K'(ro,— yn) = 
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K(xQ, + yo)i the integration under consideration can be 
expressed as 

/J A'(x(i,|j/ol) ^^ 2 / 0 = 2/ K(io,yo) dyo (B9) 

where, on the right, the absolute-value signs on ijq can be 
dropped. 


After performing an integration by parts by letting 


and 


7 o ^2/0. o 1 

dv~2 — j —7 17=— 2 — j 

I Vo I yo 


(BIO) 


iMfcyo 

ifcyo , 1 


imyo 


^ A/rp-j- \ 


I («-.vV«’+w) 


p » 

jMlfi 




Uc ^ (x-AfV^*+W) 


Ik n 
Mjo 


or 


rf«=|[- 




ikMxoVtt . 

2 ' « 2 I a2^, 2 » 


ikyo 


-f 


d\ (Bll) 


(xp® + /? ^ W Va^o^-f fi^yo‘ 


^ -^V+i^ip) _j_ 


J ’® fi—iki/ttT Pi 

^™dr+Ft/„ 

w/^Vi+T 

'I 








(B12) 


there is obtained for uv 


^ \ yoL /8 

( -^ — +iy 

Vvv+^v Mr 


lAfirys 




« (to— 3/ -^zy+0'yi^) 


+ 


i-v f" vi +T’' 

J.\///5 dijo 


d\ 


T'l 

Ju^ = oJ 

(B13) 


This expression vanishes at its upper limit yo= °° and is 
singular at its lower limit i/o=0. However, by not making 
2^0 in the derivation of K{zQ,yQ) until after tliis stage is 
reached, this singular value is canceled by other terms tliat 
have otherwise been dropped. Thus, the expression (B13) 
may be considered to be zero, which is the value of its finite 
part. The integration under consideration is then reduced 


to the value of — v du which 




I rkMxp 

+/3W"^^3-oH0V 


v*o'+/3W J .t AW ^/r+ 


ip 


, r»e^ 
Jo -y! 


-f: 


(x-A/Vv+W) 


r+/3V 


dX ydyo 


/{HI 


ik 


\ ^0 


.) 






A///»Vl + P 


rfT + F 


^X 


0 


dX 


dyo 

(B14) 


The terms of this expression are treated separately in the 
next three equations: 


First (see ref. 22, p. 180) 
^ikM . 


“X:e 


ik 


^0 J 


-gi Cxa— A/ Vxo^+^yoO j 2ik 

^ dyo^'-^ e 




ikM , , , 

— ^ koi cosh ( 


do 


H,«> (^^^')] (B 15) 


second 


and third (see ref. 22, p. 180) 


^eo ^oa j p a> 


=-2ik 


dr 


./a/ZA T -y/l -|" P 


'-2ifclog^ 


(BIG) 


"'X. ■'■'"X' 'v>?ws* '‘'‘’“X‘‘”’*X vv+»v ''‘'•- tX.'''‘'‘'X. 






iriF r^o ~ 

-- rJ. " ■ 


= //.‘"(^/|X|) 


rfx 


(B17) 
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Substituting the results in ec] nations (Bl5) to (Bl7) into 
equation (Bl4) gives 



This result is a form of the expression for the kernel function 
of Possio’s integral equation relating pressure and downwash 
for a two-dimensional oscillating wing in subsonic compres- 
sible flow. It checks the results given, for example, in 
reference 27. 

Reduction of kernel for M=l. — The kernel function for 
d/=l ma}'' be written as fsec cq. (47a)) 



(B19) 


The second integral appearing in this equation can be shovTi 
to cancel several of the terms so that the kernel becomes 


K(ro,yo)M. 


r 2 'Yj 
—/-■2 V" 

LA 

L2/0 

91 "/UA 

2/0 Vo J 


(B20) 




(B22) 


Finally, the kernel for the sonic case in two-dimensional flow 
may be written as 


l£^K{xo,yo)n,.^dy„=^^l^ 






(B23) 


It ma}' be noted that the integrals in this equation are readily 
expressible in terms of Fresnel integrals 


and 


C{x)= f cos ~ fdt 

4/0 ^ 

sin ^fdt 


Reduction of kernel for M=0 . — For M=0 it is convenient 
to modify the kernel function before integrating with respect 
to yo. For this purpose use is made of the relation (see cq. 
(B7)): 


so that the kernel for the sonic case in two-dimensional flow 
may be written as 




I K(xo,yo)M^idyo- 




f XK \ 


(B21) 


-A K,ik\yo\)-^^ [/.(iiyo|)-i,(i-|yol)] 

=1^ f ” <yIWe-‘>^^dx 

yo Jo 

and the relation 


(B24) 


Integrating equation (B21) b}^ parts with respect to yo, re- 
taining only finite parts of the integrated results, and making 
use of the relation 



—o f r ■y,^yo^-\-\^e~^^^d\ (B25) 

Vo Jo Vo J-^Q 

With these relations the expression for /iC(a:o,|/o)Af-o, equation 
(53) , can be vTitten as 


K(xo,yo)^f^o■ 



. I ik\Xo‘+yo^ 

a' yi?+^ ~ 




(B2C) 
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But 


F 


'.r J -j,. ^ ?/o' yo* j -/o ^ \2 




_ j*(r+?/o^ 

2 ' ““T „ 




til ore fore, 


^(•^o<?/o).u = 0- 


-/A-z, r 

"7^J- 




r/X 


2/0' 


(B28) 


Tntogratiiig with respect to ?/o gives 
r°° ^ 2 r® r® ^-i^x 

Jo J-.0 

“7" J-./ Jo (yo‘+xr'^ 

9 r” ^-ixx 

= _:^e-xor„ « _,;x 

t J A 

(B29) 

Integrating each integral in equation (B20) aiul retaining 
onl}" finite parts yields 

e ~ 


J « 2 / T" 

K(x„,y„).i[^odyo=—j e 

<‘r-) " 


JX- 


yo^V+2/o’ 


-•Ou 


JX 


(B27) 




il- 


/■ 

47t/' 


*+‘J,.'-x“''0 


cos ^*X 


f(-2i+X-' 

j[si(W+^]}) (B30) 


where Ci(/.*wro) xSi(^'j-o) denote, respectiveh', the ^‘cosine- 
integral” and “sine-integral” functions defined as follows: 


Ci w=-j^” dl 


V 7T sin / 

X.W-2-J, — 


(It 


The results in the braces of equation (B30) cheek with results 
given for this case in reference 14, 


APPENDIX C 


SOME REMARKS ON EVALUATION OF THE KERNEL 
FUNCTION 

Exa(“t expressions for the kernel function J5 l(jo,2/o) 
given in equation (20) for 0^d/<Cl and in equation (47) 
for i\/=i. Corresponding approximate forms are given in 
equations (54) and (56). 

Equations (20) and (47) are valid for any set of values of 
d/, A", Joj ijo. To calculate the value of tlie kernel from 
these equations, it is necessar^^ to evaluate niimcricall}" 
the integrals which appear. Values of the other terms can 
be obtained by making use of existing tables. Extensive 
tables of the Bessel functions ivi and 7i may be found in 
reference 28 and a table of the Struve function Li with second 
and fourth differences for interpolation purposes may be 
found in reference 29. Sample values of the kernel are 
given in table I. 

For certain ranges of values of d/, 7, a*o> 2/o» as indicated 
by equations (55) and (57), the kernel can be evaluated 
hy making use of the power series expansions given b^’ 
equation (54) for 0^d/<l and equation (56) for 47—1, 

The various expressions for K(jq,ijq) become singular when 
^^= 7 /— In order to be able to evaluate the kernel 
in such circumstances, it has been separated into two parts 
as shown in equation (30). One of these is denoted b\" 
K(xo,yo) — K' (xo^yQ) and is not singular; the other is denoted 


b\" K'^Xo^yo) and contains all the singularities. Obtaining 
the value of (K—K') from the form of the expression given 
in equation (30), however, may be troublesome. This 
particular value for yo^O, Jo> 0 can he obtained from the 
following limiting form : 


lim [K(xQ,yo) ■ 

Vd-»0 

-K'(xo,yo)] = 

j 

^~P~\ 


a-(l+2\/)\ 

2xo / 




ik.k- 

r2+M 

Li+i»/ 

-27— log ( 

kxo \ 

i+Mr 


Ci ( 

( kxo \ 

<1+Mj 

i+{ Si ( 

kXu \ TTl 

1+4/y 2 

]} 


where y denotes Euler’s constant (7 = 0.577216) and Ci and 
Si denote cosine-integral and sine-integral functions, re- 
spectively. (Tliese functions are tahulated in reference 30), 
For 4/=l, this expression reduces to 

lim [K(:ro,yo)~A’'(.ro,yo)]-%-4 

•' Xu 

I [|-2T-log(^^”)+Ci (-|“)+i Si (C2) 
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The kernel function is not sing:nlar for For ?/y=0 

and Jq< 0 it may he written for d/<l as 

i/c-K) f L L--^0 J 

l[<-'KS/)-'KSr)+?]} p‘' 

The expression for K—K' for .r„ CO, ?/d = 0 may also he useful. 
It is 

lim [Kic-To;y^-K’{-x,^i^^''p-{ 

yc-0 ( [^0\ 

\2xo^ 2\xo\ J 

5r^r/-'"KS7)+<^<S/)- 

For K{~xo,yo)=K{—xo,yo)~K'{—jo,yo) = 0. 

Some results of evaluating the kernel and its nonsingular 
part are given as examples in table I. (In order to obtain 
these results the required integrations were performed 
numerically by manual computing methods.) 

APPENDIX D 

ALTERNATE FORM OF EQUATION (20) 

Subsequent to the derivation of equation (20) as given in 
the textj it was found that the two integrals involved in this 
equation can be combined in a manner that leads to a more 
concise and, for many purposes, a more convenient form of 
expression for the kernel function. The purpose of this 
appendix is to derive this alternate form. 

Consider first the integral 

Gi=TTfn — (Dl) 

A/(i-yo)Uo 


and make the substitution 


i (X-A/vV+^^a-yo)0 = -^i?/olr 




X=/-|2/„|(AATT^-r) 

This suhstitiition gives for Gi 


(D2) 

(D3) 






Mk\y,\hri0 hl + r^ 7 

1 r % , 

liVoiJ.A///? 






k'l-r r- 


Consider now the integral 

C3//I8 ^ 

G,=j 


(D4) 


(Do) 


and integrate by parts l)y letting 


?/ = ^ ! — T- 

<lc=f~‘^' ^ ^ (It 


so that 


rhi- 




v= f -, 
I //m1 


This integration gives for Gy. 






^ I h ^ - p-iKif 'r,]r 

^lh\ Jo ^1-i-r' 

Subtracting G. from Gy (eqs. (D4) and (DO)) gives 

f7,-f7,===,7-j---7^ g“ 0 ^ ® "+ 

^‘i?/o 0k\yo\ Mk-y^ 


(DO) 


1 ^ 


Alk'yQ 




+ 


i T 

. I , e’ 

^ i?/0i */ 0 1 -f~ T^ 


i '^Vzn-'+ZJ^KoO J . 

^i?/o| Jm/0 




7-i2/„] ^mkyof '^m-yor 

7-i (zo-MVz«-+W) t 

— ^_z=. (D7) 

^l2/o| .io A i 7 


4- p 

^ij/ol Jo 


Substituting this result into equation (20) of the text gives 
for K{x„,yo) 

A>„, [7,^ 

i kxo ~ [kzo-^r^{kI,y--\~^^(.ki/o) 


!/o I {kyo)^\(kxof (kyof 
i (rc-.Arv^c-'+w) r r , 




1 + 7 


•} 


+ 

(D8) 


The integral in this equation is in general more amenable to 
numerical evaluation tlian either of the two integrals appear- 
ing in equation (20). Furthermore, with this expression, it 
is not necessary to consider tlie incompressible case as a 
special case, since no trouble arises in setting d/=0. Simi- 
larly, for the sonic case no trouble arises and this expression 
gives for Jq>0: 

A( J’o, 7 /o)a/ - 1 = J2 = yy ! ^ ~ 


1 J f; 

iiyolJo vi+7 J 


+ 


(D9) 
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TABLE I.— ALVLUES OF THE KERNEL AND ITS XOXSIXGULAR PART AT J/-0.7 


1 

tq 

Z/« 

‘ 

^ (-TojI/o) 

K{xo,yQ)~K'{xo,7jQ) 

0 

0. 125 

0- 1 

-63. 827569+ I. 112400z 

0. 1 44529 007824/ l 



. 3 

- 63. 801759+ 3. 290793i 

003441 069879/ 



. 5 

- 63. 5 1 3049 + 5. 408465/ 

00942.3-. 192655/ ! 



. 7 

-63. 127659+ 7. 466762/ 

0181 14-. 374807/ I 



1. 0 

-62.396691+ 10.445693/ 

035609- . 756548/ ; 

1. 5 

. 125 

. 1 

-126.263912+ 19.142811/ 

, 141754-. 028841/ I 



. 3 

-114.855158+ 55.631898/ 

-. 0313 17-. 056060/ 



. 5 

-92.964383+ 86.829346/ 

-. 1 23447-. 115703/ 



. 7 

-62. 878740+109. 92702(h: 

-. 283001 133318/ : 



1. 0 

-8. 792808+125. 223964/ 

-. 581313 + . 022309/ | 

0 

6. 0 

.1 

019271+ , 016639/ 

000039 -. 000699/ 



. 3 

+ . 007493+ . 020950/ 

+ . 007793-. 049064/ ^ 


! 

. 5 

+ . 020861+ . 001545/ 

. 03G1G5-. ] 15145/ 



. 7 

+ . 009570- . 017888/ 

. 095337-. 181254/ ' 



1. 0 

018833- , 006290/ 

, 305G27-. 239670/ | 

1. 5 

6. 0 

. 1 

-.027209+ .020038/ 

00905 -. 006215/ 



. 3 

+ .002452+ .028186/ 

005415-. 041401/ 



. 5 

+ . 021871+ . 013305/ 

007432-. 109920/ 



. 7 

+. 022588- . 008980/ 

-. 026790-. 232786/ 



1. 0 

-. 004786- . 022987/ 

-. 190134-. 523276/ 
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